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Abstract 

It  has  been  shown  that  literal  information  can 
enhance  geolocation  information  in  the  multi-target 
tracking  and  data  association  problem.  This  paper 
continues  previous  efforts  in  establishing  a  sys¬ 
tematic  approach  to  the  combination  of  both  types 
of  information  using  membership  functions  based 
upon  multiple-valued  logic.  Filters  are  established 
for  literal  and  non-numerical  attributes,  somewhat 
analogous  to  the  well-known  Kalman  filter.  The 
major  result,  however,  is  an  improvement  and  clari¬ 
fication  of  a  previous  theorem  establishing  asymp¬ 
totic  forms  for  the  posterior  possibility  distribu¬ 
tion  of  the  unknown  data  association  parameter  as 
Information  granularity  decreases  and  as  inference 
rule  structures  become  more  definitive. 


1.  Introduction 

.The  multi-target  tracking  and  data  association 
(or  as  commonly  called,  "correlation")  problem 
still  remains  the  center  of  much  activity  and 
interest.  In  the  past,  emphasis  was  placed  upon 
the  use  of  only  geolocation  data-i.e.,  information 
containing  reports  on  (usually)  two-  or  three- 
dimensional  target  positions,  together  with  possi¬ 
ble  velocities,  accelerations  and  related  equations 
of  motion  parameters  , (|1 -3] .  More  recently,  an  ef¬ 
fort  has  been  carri^' out  in  utilizing  in  a  more 
rigorous  or  sys tema tic  manner  other  types  of  inform¬ 
ation  ,  including  various  sorts  of  sensor  system 
parameters  that  could  also  be  formally  treated 
in  terms  of  equations  of  motion.  In  addition,  use 
of  attribute  information  of  radically  different 
natures  is  also  sought.  This  includes  visual  sight¬ 
ings,  classifications,  end  other.often  non-numerical 
and  linguistic-based  (or  narrative), information 
which  cannot  be  treated  from  traditional  statistical 
analysis.  This  work  is  presented  in  [4-6],  based 
upon,  and  related  to, in  general,  previous  and  on¬ 
going  research  in  possibility  and  fuzzy  set  theory 
and  its  real  tionships  to  classical  >.rou<ibility 
theory  and  multi-valued  logic  and  set  theory  [7-9]. 

In  this  paper  we  continue  to  attempt  to  estab¬ 
lish  a  unified  approach  to  the  integration  of  both 
types  of  information  for  the  multi-target  tracking 
and  data  associat+on  problem.  However,  since  the 
results  are  carried  out  on  a  general  level,  with 
suitable  modifications,  applications  to  medical 
diagnosis,  fault  determination,  and  other  problems 
involving  non-sequential  knowledge-based  systems 
can  be  established  [9]. 


As  stated  before,  this  work  still  unfortunate¬ 
ly  requires  both  a  sense  of  "art"  as  well  as  sci¬ 
ence  in  achieving  its  goal. 


2.  Basic  Problem 

The  basic  problem  can  be  stated  as  follows; 

Observed  or  reported  data  arri ves, usual  1 y  in¬ 
dexed  by  time  t^^  ,  attribute  type  k,  and  sensor 
source  u.  For  purposes  of  simplicity,  we  will  omit 
the  last  index  and  assume  the  multiple  sensor 
source  problem  is  resolved  and  plays  no  role  in  the 
analysis  here.  Attributes  A,,  k  =  l,2,..,H",  coo  be 
either  objective  or  statistical  in  nature  or  they 
can  he  subjective  or  linguistic  in  nature.  Denoting 
the  natural  domain  of  values  (which  can  be  numbers, 
vectors  of  numbers,  linguistic  labels,  etc.)  of 
each  attribute  A.  as  D,-  dofn(A,),  statistical  att¬ 
ributes  typically  have^their  domains  D.CIR'^k  with 
error  distributions  in  the  form  of  classical  prob¬ 
ability  functions,  typically  being  discretized 
versions  of  Gaussian  or  mixtures  of  Gaussian  dis¬ 
tributions.  Thus,  in  general,  unless  representing 
dirac  or  mass-point  distributions,  the  error  prob¬ 
ability  functions  are  not  normabl e , 1 .e . ,  they  have 
maximal  values  less  than  unity  (due  to  the  obvious 
constraint  of  probabilities  adding  up  to  one).  On 
the  other  hand,  subjective  attributes  Aj^  typically 
have  D.  being  some  finite  set  cf  labels  with  no 
natural  spacial  ordering  present,  unless  D.  repre¬ 
sents  a  sst  of  n|^-tuples  of  measurement  quantities, 
where  for  example,  A-=  "large"  could  have  D^  con¬ 
sisting  of  pairs  (a,D),  with  a  being  weight  in  lbs 
and  b  length  in  meters.  The  error  distributions 
corresponding  to  subjective  attributes  are  often 
in  the  form  of  normable  possibility  functions  -  but 
not  probobility  functions  -  obta ined  from  a  panel  of 
experts,  in  contrast  to  the  geoteation-physical 
derivation  of  statistical  attribute  error  distri¬ 
butions.  For  motivation  for  such  subjective  distri¬ 
butions,  representing  overlapping  and/or  vague 
compound  events  -  in  contradistinction  to  statisti¬ 
cal  uis*'r’butions  representing  disjoint  exhaustive 
events-  see  again  [4].  Each  error  ..liti  ibution 
represents  the  conditional  possibility  (possibili¬ 
ties  include  probabilities  as  special  cases)  that 
a  particular  attrih"«-p  Hoo^iair,  ,-3''uc  is  3,.tj3llv 
present,  given  an  observed/reported  domain  value. 
Whether  statistical  or  subjective,  it  is  assumed 
that  at  least  theoretically  each  such  distribution 
is  obtainable. 
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In  addition  to  error  distributions  and  observ¬ 
ed  data,  a  collection  of  inference  rules  is  assumed 
to  be  present.  Each  such  rule  R  consists  of  an 
antecedent  part  ant  and  a  consequent  part  conseq 
connected  by  :n  impYication  operation,  ant  con¬ 
sists  of  the  conjunction  (more  generally,  sXme  com¬ 
bination  of  conjunctions,  disjunctions,  and/or  ne¬ 
gations  )  only  of  modified  matching  tables  for 
each  attribute  from  some  subset  of  all  the  attri¬ 
butes.  Each  modification  mod^  of  basic  matching 
table  M.  for  A.  is  assumed  ’  to  be  in  the  form 
of  an  intensification  or  strengthening  or  an  exten- 
sification  or  weakening.  More  general ly ,modi fica- 
tions  can  involve  positive  and/or  negative  forms. 
For  simplicity,  we  assume  throughout  here  that  only 
positive  forms  are  present.  (See  [10]  for  more 
background  on  modifiers  and  hedges  for  attributes.) 
In  ordinary  English,  intensifications  can  be  rep¬ 
resented  in  a  simple  order  as  e.g.  ..."slightly 
very" , . . ."very" .... , "very, very" , . .  ."extremely"  .... 
Extensi f ications  can  be  expressed  as  e.g.  ..."very 
very  little  of".. .."very  little  of". ..."more  or 
less  of"  ...  ."moderate  amount  of"...  The  conse¬ 
quence  consists  of  a  single  modification  mod 
(again,  either  a  positive  intensification  or'^ex- 
tensi  ficatlon)  of  the  basic  data  association  level. 


Then,  as  data  arrives,  track  histories  are 
built  up  (with  some  eventually  discarded)  sequent¬ 
ially  in  time,  consisting  of  those  data  reports 
which  are  adjudged  oS  belonging  to  the  same  target 
source.  The  decision  on  associating  together  or 
not  a  given  pair  -  a  previously  established  or 
tentative  track  history  with  a  new  report-,  is 
based  upon  the  application  of  the  set  of  inference 
rules  and  error  distributions  upon  the  two  data 
vectors,  where  it  is  understoocLtbat  the  previousl.v 
established  track  history  i  is  suitably  updated 
for  comparisons  with  the  new  report  j.  Of  course, 
for  a  given  track  history  the  updating  of  the  data 
history  for  a  statistical  attribute,  under  Hnear- 
Markov  or  Gaussian  assumptions  .can  be  inade  using 
the  now  well-established  Kalman  filter  [11].  On  the 
other  hand,  the  very  different  structured  subject¬ 
ive  attributes  require  another  approach.  For  ex¬ 
ample,  A.,  might  be  "color  of  flag"  with  0.  =  {red, 
red-orange, orange  .... ,  blue, ... .black  striped,..), 
or  A.,  might  be  "class"  with  0.  =  {  C,  ,C_,...). 

How  00  we  "filter"or  predict  among  such  Values? 

Thus  an  attribute  filter  is  sought  for  these  cases. 
But  after  a  moment's  thought,  the  answer  to  this 
problem  is  relatively  simple:  obtain  the  possibil¬ 
ity  analogues  for  the  p'obability  function  situa¬ 
tion.  This  is  spelled  out  in  section  4.  (  For 

a  schematic  view  of  one  cycle  of  the  data  associa¬ 
tion  process  from  one  sampling  time  to  the  next, 
see  Figure  1  at  the  end  of  section  4.) 


3.  Definitions  S  Notation 


T*-"  -otaiion  used  here  is  somewhat  different 
from  the  previous  and  related  paper  [6],  due  to 
the  ongoing  effort  to  emn’oy  the  sinl^llcJt,  yet  nictl 
accurate  notation. 


o(i.J)  i  fo(i.J)\ 

i  ■  ^®v  v=l,...N 


€  dom(ci’’J^)  =  [0,1]  is  unknowri  parameter 


representing  data  association  level  between  trad 
histories  i  and  j  ;  the  attribute  represer" 

"correlation"  between  i  and  j. 

Aj,A^,..,A^,  is  the  set  of  attributes  which 
are  statistical  ; 

+] +2  •  ■  -  5et  of  subjective 

attri  butes. 


To  indicate  dependency  and  chat.ge  due  to  samp 


inq  time  t  ,  where 
m 

V  '3^ 


dom(  A,  ^  ’  ^ )  =  D,  ^  ^ 
k ,m  k  ,m 

similarly  replacing  rj|^; 

iX(D  . 

m  ,  A  K  ,m  > 
k  =  1 


'k  ,m 


k  ,m 


ii) 
'k  ,m 


"^k  ,m 


updated/smoothed'observed  values  for  A,  , respect 

1ve1y.  ^ 

2  v'i.j)  A  f?  ('1  7  (j)i 
H.m  ■  '^k,m  •'^k.m  '  ’ 

\  >m<j„  .  for  any  m^  =  0.1,?,., 

(i.mr,)  i  (?'i  ))  °  , 

0  m  m<m.  » 


Z 


z  i  (1  (’h 

m  ^^k.m  ''k  =  l  ,..,M"  ’ 


etc . 


Similar  notation  holds  for  the  true  possible 
attri bute  values . 

We  use  primeC )  notation  to  indicate  the  stat 
istical  part  »nd  double  prime  notation  (")  to  indi¬ 
cate  the  subjective  part: 

7.(1)  i  ,7  {i)i 

m  '■'"k.m  'k  =  l  ,.  .,M' 

7..(1)  i  ri  (i)\ 

m  ,m  'k-M'+l  , . .  ,M" 


p(i )  ^  D'  ^ 

m  m  m 


,  etc. 


We  use  ♦  to  indicate  possibility  function; 

♦  (■I")  to  indicate  conditional  or  dependent  poss 
ibility  function. 

„  (i: 


k  ,m 


is  the  index  denoting  the  conditional 
(i)  .  Thus 


error  possibility  function  for  A 


and 


♦  (■l-dPp])  -’[0,1] 

K  jTH  k  I  fn  k 

♦  (Z  (i);pj’^) 

k  ,m  '  k  ,m  k  ,m  ' 

=  possibility  Z,(1)  is  true  value  for  Af’2 
k.m  k.m 

g1  ”  'is  nnditrd  uuserved  "aiue. 
n.  ,ni 

Also, 

p  ( i  •mp )  4fp(i)\  _  /p'(i  iiHo  ^  p  ”  ^  f  'ifio  ^  i 

®k  ‘^^k.m'msm^  '  '  k  '^k  ' 


p(0  =  (p'(i)  p..(i) 

m  m  ’  m 


^"^^k.m^d  , . .  ,M"  • 
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Four  special  possibility  functions  are  indi¬ 
cated  by: 

0  :[0.1]  [0,1], a  decreasing  or  at  least 

nt 

non-increasing  function  with  boundary  conditions 
coinciding  with  usual  two-valued  logic; 

l=fnt'05  ^  0=*nt''’  • 

representing  negation  or  general  set  or  attribute 
complement  ; 

♦  [0,1]  -^[0,1],  a  nondecreasing  function  in 

o 

each  of  its  arguments,  which  is  conti  nuous  min , 
associative,  symmetric,  and  possesses  the  boundary 
conditions  reducing  to  the  two-valued  logic  ones: 

0=it,  (x  ,0)  =  *,  (0,x)  ;  x  =  iti,  (x  ,1  )  =  ♦,(!  ,x) , 
for  all  X  [0,1]. This  operator  may  be  unambig¬ 
uously  extended  to  an  arbitrary  finite  nur.ber  of 
arguments  and  represents  conjunction  or  general 
set  or  attribute  intersection; 

2 

[0,1]  -^[0,1],  a  nondecreasing  function 

in  each  of  its  arguments,  which  is  continuous, > 
max,  associative,  symmetric,  and  possesses  the 
boundary  conditions  reducing  to  the  two-valued 
logic  ones: 

for  all  X  E  [0,1].  This  operator  may  be  unambig¬ 
uously  extended  to  an  arbitrary  finite  number  of 
arguments  and  represents  disjunction  or  general 
set  or  attribute  union; 

2 

♦^:[0,l]  -►[0,1],  a  non-increasing  function  in 

its  first  argument  and  a  nondecreasing  function  in 
its  second  argument  with  the  two-valued  logic 
boundary  conditions  —  - 

)=<i^{0,x)  =  ;  0  =  <>^(1  ,0), 

for  all  X  E  [0,1].  This  operator  represents  logic¬ 
al  implication  :  "if(-)  then  (••)"• 


is  also  called  in  the  literature  a  t-norm 

and  (f  is  called  a  t-conorm.  See  [12]  or  [9], 
Chapt?l"  2.3.6  for  more  details  on  these  function 
classes . 

Next, let  denote  the  inference  rule, 

with  associated  possibility  function 

*(-|--;R^):[0.l]  -  [0,1], 

somewhat  abusing  notation,  where  there  is  an 
associated  index  set  {l  ,2, . .  ,M" ) ,  where  as 
usual  , 


so  that 


•■nt  <'rl  "”>■  ‘s/'nod  <*K 


♦  (e  4-  =  i()  .  (ir-  (b  ))  . 

^conseq^  v'  mod,  ^Co  v" 

Denote  also  the  index  set 


Then  we  write  for  the  overall  error  distribu¬ 
tion  : 

(7(i.j)|^(i.j)  p(>.j)^ 

))^ 

4(Z'[''j^l5'(ij);p'' 

m  ‘  m  m 

i  *  (  a(7^'  •J'n- 

'/r  '%,m  ’^,m 

(k-1 


k,m’  k  ,m 

?" 

'm  •  m 


with  p{  Z"( '  '  3  ]  I  ’-i  ^ )  defined  similarly 

m  •*' 

in  terms  of 


For  the  overall  inference  rule  effect,  we  have: 


rv=i . . .  ,N^ 


R  1). 


(?) 


Other  notation  and  definitions  will  be  intro¬ 
duced  as  needed.  From  now  on,  only  the  most 
important  changes  in  t-norms  and  t-conorms  will  be 
specially  indicated;  otherwise,  the  generic  notation 

,  4  wi 1 1  be  used  . 

8.  or 


■4.  Basic  Analysis 

The  following  simple  theorem  is  proven  easily 
using  the  basic  properties  of  conditional  possibil¬ 
ity  functions  as  given  in  [6], [7],  or  [9]  and  car. 
serve  as  the  desired  attribute  filter.  Note  also 
that  in  general  the  distinction(omitting  P,(i.>'’)) 


Theorem  1 

Suppose  (omitting  the  andii'^  1  for  any 

k,  k=M'+l ,.. ,M", 

are  known  functions  with  4(7^, 

normable  with  respect  4(7,,  and  without 


loss  of  generality  suppose  that  4(7  |Z. 

normable.  Then-  ’ 


k  ,m '  ‘'k 
e  t 

(*,(P(7 


(m) 


|0(m-l ) 
k ,m ‘ ^k 


)  is 


k  ,m 


'^k,m-l 

;(m) 


(•j<*lZkjjl!"‘'''«<VlZk*») 


is  implicitly  obtainable  from 


(m-l  ) 


■•.'•<Zk,.ll£"'’>.‘'lk,JZk,, 


)) 

)). 
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Next,  ^or  purpose  of  completeness,'  t^e  full 
fACT  algorithm  [A], [6], [9]  (Chapter  9)  will  be 
presented,  for  justification  of  the  results  see  the 
above  references. 


The  fi( 
is  4 ( 6  I d i a q , 


jr  possibility  function 

•J  ^  )  obtainable  implicit- 


ly  from  the  equation  ,  again  using  simple  notation, 

tjiggC-  I  -r,, 

(t. 


^  ;R,P)=4,('4(e|diag,?,l’J  ^  ;R,P), 


*.  (d'''ql“ 

'univ  ' ^ ^  >^"diag 

(Os.esl) 


(diag)),(3) 
(elZ^J^  ;R,P)),  (4) 


■.K.PUs 

v  =  l'^.  .  .N 


^;R,P)=4^^(4(e.zi’^l^'J^R.P)).  (6) 
(l  (ij)cD^'J^) 

\  rn  m  ' 

4(e,Z^^J>|!_[^’J>;R.P)=4j(4{6|Z^'j\z^’j^R.P),  (7) 


where 


4(e|Z^''^^'i^’J^R,P)  =  4(e|Z^'-^’;R)  .  (8) 

*(z'’JVzi'j^R.P)  =  ♦(zi'^’|zi'j^P)  (9) 

are  obtainable  as  in  section  3. 

Furthermore,  in  general,  4  above  breaks 
up  into  statistical  and  subje??ive  components: 

4(6|Z^'j^',R,P)-4^^„  (4„..(G(6,zi’j  ••,«,?))) 


where 


G(0,Z^'^'^;RP)-  —  - 

\(*H(e,zti.M:z;(ij);R.p'‘).4(z^(ij)|?;,(ij);PH 

4j(4(e|Z^^J^R),  4(Z;(''J)|?^^'J’;P")),  (12) 

where  4^^,,  may  be  chosen  as  max. 

Since  in  general  ’5  a  nondecreasing 

function  in  e  with  value  1  at  6=1,  and  if  all  modi¬ 
fiers  are  chosen  in  the  positive  sense  as  discussed 


in  section  2  so  that  each  4 


mod,  Vod 

k  ,v  v 


are  nondecreasing  functions  of  x  with  unity  values 
at  x=l ,  such  as  is  the  case  with  exponentials, then 
it  follows  from  the  property  of  4^  that  the  final 
posterior  possibility  function  is  formally  like 
a  distribution  function  except  for  not  necessarily 
being  0  at  6=0.  Hence  a  reasonable  measure  of  cen¬ 
tral  tendency  for  0  so  described  is  ,  using  formal¬ 
ly  expectation  notation^ 

E(e|diag,?^^’J^;R,P)4  je.d4(e|diag,Z^’-^^R,P) 

6=0 

=  1-  J  4(e|diag,zJ,’^^R,P)-de. 

,  0*0  HTi 

(See  also  [6]"for  further  discussion.)  '  ^ 


A  simplified  outline  of  a  data  association 
procedure  utilizing  the  PACT  algorithm  as  describ¬ 
ed  above  is  given  in  Figure  1. 


Figure  1.  Flow  Chart  for  Data  Association 
From  One  Data  Sampling  Time  to 
the  Next 


5.  Some  Asymptotic  Results 

Given  the  above  scheme  for  obtaining  the  final 
posterior  possibility  function  for  the  level  of 
data  association  between  two  track  histories,  the 
following  natural  question  arises:  Can  some  analogue 
with  the  classical  case  of  statistical  consistency 
be  established  here?  In  [6]  it  was  shown  that  as 
the  fineness  of  the  domains  of  statistical  attri¬ 
butes  increased-  and  dually,  the  granularity  de¬ 
creased  -  under  certain  reasonable  conditions, the 
posterior  function  for  data  association  converged 
to  a  computable  function  of  a  statistical  expecta¬ 
tion.  In  turn,  the  behavior  of  this  limiting  ex¬ 
pression  as  inference  rules  and  error  distributions 
became  more  precise  was  considered,  leading  to, 
under  perhaps  too  stringent  conditions,  a  form  of 
consistency.  In  order  to  clarify  this  further,  a 
similar  result,  under  weakened  conditions  and  with 
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more  specific  applications,  will  he  presented. 

Si'ppose  the  following  additional  assumptions 
hold  relativeto  the  basic  situation  used  to  deter¬ 
mine  the  posterior  function  of  6; 

(a)  for  all  y  t  .  f[^],(-|y):  R"''  -  r"'  is  a 

bounded  continuous  conditional  probability  density 
function,  for  k  =  l  , .  .  ,M' ,£=1 , j  ,  for  some  chosen 
track  hi  stories  i  , j  . 

(b)  For  each  integer  p>l  ,  let  be  the  p*"^ 

discretization  and  truncation  of  D-kS-lwith  a,  ‘f/ 

,  k,m  k,m,p 

denoting  the  mesh  of  D,  1'^'  (i.e.,  maximal  length 

k,m,p 

of  any  rectangle  formed  within  the  truncation 

area),  so  that, in  any  sense, lim  D,  vf)  =  R^k 
’  ’  k ,m ,  p 


assumptions  id)-if'  hcl.!.  Thpr  ng  ;  ,,-max. 


(f) 


p-^-foo 


lim  A  =  0  ;  in  addition,  replace  each 

k,m,p 

p-va^oo 

)  by  the  approximating  probability 
function  (I'l) 

*  (Z  ^^'l=f 

'p  ^k,m  '^k,m’  k,m  'k,m  '^k,m'  k,m^lynp 

and  denote  accordingly  all  computations  involving 
this  replacement,  for  k  =  l,..,M'  ,  £=i,j,with  p. 

(c)  '5  analytic  about  x^  ^=0, 

lt  =  l  ..  ..M'  ,f  =  i  ,j  . 

(<i)  *or' Archimedean  t-conorm  with  genera¬ 
tor  function  h:  [0,1]  -«  IR^U  {+«>)-which  is  continu¬ 
ous,  non-increasing,  with  h(l)=0  and  h(0)  < 
and  is  such  that  -  - 

♦nr'  (n'in(£h(l-x  ),h(0))), 

for  all  X  c  [0,1],  q  =  l,2,..  .  (See  ,e .  g.  ,[1 2]  and 
eg. (23).)" 

(e)  Referring  to  (d),  h(l-4  (x,y))  is  analytic 
about  x=y=0  in  [0,1  ]^. 

(f)  ifie  I  Z' ’ ^  ;R)  is  a  continuous  function  in 

Z^{’«J1  allowed  to  be  arbitrary  i  n  XF^*^)  . 

(This  is  guaranteed,  if  is  continuous  in  its 
first  argument  and  is  continuous  over  [0,1] 

for  k=l,..,H'  ,  v  =  l,..,fe:'') 

Next,  define  , 


k  =  l....H'  f=T.j  / 

’'k,r°  • 

k=1 

and  (16) 

u(x)  ?  -(d  b(s)/ds)j,,,  •(34.j(x,y)/3y)^^Q  , 

for  all  X  t  [0,1]. 

Theorem  2 

Assume  the  basic  situation  holds  as  establish¬ 
ed  in  the  previous  sections.  Suppose  also  that 


(Ifi) 


(^(f.z"^' ’j^))  (iR) 

-  m  , 


^  lim(  ♦n{o|Z^’  '^^R,^)  ) 
p  —  tn 

P  ^4  00 

=  1 -h'^ (min(p{t ] ,h(0) ) )  , 
where 

o(e)  =  max 

(Z^(i.j),  o'Ki.j)) 

a(e,z;(’’'’’)  =  r.Ey(v(c^(e,V,Z^<''-J^)))  ,  (?0) 

where  fy{-)  denotes  ordinary  statistical  expecta¬ 
tion  with  respect  to  random  vector  V  (replacing 
the  argument  non-random  Z^fi.j))  which  has  prob¬ 
ability  density  function  f  where  at  any  value 
Z'(i»j)  for  V,f  has  the  conditional  form 

^  =  1  ....M' A 

V  f=1,j  ) 


(21  ) 


(Proof : 


Use  the  canonical  expansion  for  as  in  (15) 

where  q  is  replaced  by  a  summation  over 
and  X  is  replaced  by  G  (e ,zi^' ’■^^ ) .  In  turn,  expand 


expand  h(l -0^(0  ,Z|[' ’■^  ^ ) )  in  terms  of  variable 
♦  (Z'(i*J)|Z  '(i'j))  around  0,  obtaining 

h(i-Gp(e,zi’-jh)=v(*,(^zl’>jh).*^(z^f'jfc«i) 


p  —  m 


it) 


Also  expanding  for  small, 

-TTC., ',(:„)  • 


where 


m,p 


("’bi:]'-) 


The  result  then  follows  from  the  definition  of  an 
integral .  ) 

Remark . 

An  important  family  of  t-norms  and  t-conorms 
is  due  to  Frank  [13].  (See  also  [12].)  These  satis¬ 
fy  the  basic  modular  relation 

♦pf.(x.y)  =  X  +  y  -  ♦j(x,y)  ,  (ZZ) 

for  all  x,y  c  [0,1],  Frank  has  shown  that  this 
family  can  be  characterized  by  the  Archimedean 
class  -  i.e.,  all  ♦  such  that  for  all  x  c 
(0.1).  * 

♦  j(x,x)<x  and  itijj^(x  ,x )>x  , 
which  is  also  DeMorgan,  i.e. 


(23) 
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v*g^(x.y)  =  1  -  . 


(24) 


<ti 


for  all  x,y  c  [0,1],  and  a'  well  by  the  class  of 
all  ordinal  sums  (  typesof  affine  transfoems 
involving  the  block  diagonal  parts  of  [0,1]^  -see 
[0].[1 2]  ,[1 3])  of  subsets  of  the  Archimedean  class. 
Indeed,  the  Archimedean  class  of  Frank  contains 
many  common  t-no,rms  and  t-conorms  and  can  be  con¬ 
veniently  parameterized,  using  parameter  s  ; 

wi  th  <t  (x,  ,..,x  )  determined  from  the  DeMorgan 
or  1  q 

or  modular  relations,  for  all  x.|,..,x  c  [0,1]  , 
with  generator  function  given  as 

h^(x)  =  -  log^s^-l  )/(s-I  ) ) ;  X  c  [0,1],  (26) 

for  all  0<s<+“  ,  where  the  special  cases  are  ; 
s=0  (non-Archimedean )  : 


h,0  =  min(x^,..,Xq) 


^or,0^>'l'---’'q^  ^ 

6o{x)  =io  , 


max(x, , . .  ,x  ) 
if  x=0  ^ 


(27) 


s  =  l 


if  0<x<l 

■  X 


s- +« 


9j,i  v*l  ..  ..^q/  -q 

^  prob5um(x,  , . .  ,x  ) 

^  1  -  n  (lix  ) 

1  CO 

i 

h.|(x)  =  -  log(x)  ;  all  x  c  [0,1] 
min(x^  +  ..^x^.l). 

*or  ,+”^^1 
h. 


(28) 


Xq)=max(x^+- •+x^-(q-l),0)  (29) 
(x)  =  1-x  ;  all  X  £  [0,1], 


It  follows  that  if  all  t-norms  and  t-conorms 
used  in  the  computations  for  the  posterior  function 
of  6  are  culled  from  Frank's  family  and  the  assump¬ 
tions  (a),(b),(f)  hold,  then  Theorem  2  is  valid, 
and  the  key  computations  for  <  and  v  are  : 


For  generator  h^,  of  ♦ 


(d  hj,  (x)/dx)^^.| 


For 


'S.s- 


(3  *j,3(x.y)/3y)y^0 


or', s'  ’ 

f-(log  s  ' )  •s’/(s'-l ) 
(0<s'<+<”,s'fl  ) 

(s'=l)  .  (30) 

(s’'-l)/{s-l) 
(0<s<+<»,sf  1 ) 


X 

(s  =  l) 


(31) 


For  ♦ 


S3,s' 


.s' 


2M'-1 


[((log  s')/(s'-l)) 

(0<s '<+”,s 'fl ) 

(s'=1)  ,  (32) 

Finally,  let  us  consider  briefly  the  effect 
of  the  choice  of  inference  rules  and  the  accuracy 
of  the  error  distributions  upon  the  asymptotic 
posterior  function  for  6  as  given  in  Theorem  2. 

So  far,  matching  tables  have  not  been  discussed. 

In  general,  although  the  basic  structures  of  the 
Inference  rules  do  not  depend  upon  time  and  the 


particular  pair  rf  trarl  histones  in  question,  tip 
matchino  table  docs.  Tor  subjective  attributes, 
these  could  be  chosen  as  some  type  of  s  yir.ri.e  t  r  i  ?  a - 
tion  o'"  the  correspond  i  ng  error  distributions, 
while  for  siutislical  attributes,  the  matching 
tables  can  be  naturallv  chosen  as 

,  1;)! 

k  ,m  k  ,m  x  ,m  tm  km  k  ,m  k  ,m ' 
and  F  is  the  probability  d  i  ;tri bution  function 

'^k  ( T) 

of  \  corresponding  to  the  7'  being  r.v.'s  with 

a  comnon  expectation,  hut  otherwise  described  by 
f,  yf)  as  for  V  in  Theorem  ?. 


independent  p.d . f . ' s 


Note  also  by  inspection  that  4^  (^,  •^  ^ )  is 

a  non-dg(^reasing  |‘un^tion  in  the  variables 


Kx  (  Zi ’3  M  ,..,>(  Zi,  )) ,  provided  that  we  choose, 

Im  n  ,m 

for  example,  *  .  (xl^x^x'.k  and  $  .  (x)  =  x^v 

mod  ,  mod 

V  ,k  V 

for  some  constant  positive  a^^,  ,  b^  ;  all  xc[0,l]. 

It  then  easily  follows,  that  if  all  error  distribu¬ 
tions  are  not  in  dirac-form,  but  indicate  some 
spread  or  confusion  between  observed/smoothed  data, 
then  if  we  can  legitimately  obtain  for  the  infer¬ 
ence  rules  the  consistency  structure  given  by 
all  a^i^  approaching  the  extreme  extensification  0, 
with  'all  b  approaching  the  extreme  intensi¬ 
fication  +»  \  for  0<e<l  ,  *(elz(ij';R)  approaches.. 
0  for  fixed  values  of  £  ,  and  similarly  for  *^(^^^'7 
On  the  other  hand,  for  e=l  ,  the  latter  is  maxV 
mized.  The  consequence  of  this  .applying  standard 
inequalities  to  the  expectation  Fw(-)  in  Theorem  2 
and  making  the  additional  mild  assumption  (true 
for  Frank's  family  for  s>l)  that  v(x)  is  non-de¬ 
creasing  1n.x<over  [0,1]  with  v(0)=0,v(l  )=1  ,is  that 
<i>(ejdiag,?y3^;R,P)  approaches  the  possibility 

function  (f  representing  complete  data  association 
between  i  and  j,  i.e.,,  . ^^3 


ff  o<e<i  . 


(35) 


Future  work  will  be  directec  toward  extending, 
further  quantifying,  and  establishing  empirical 
bases  for  the  above  study. 
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